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By means of full exact diagonalization, we study level statistics and the structure of the eigenvec- 
tors of one- dimensional gapless bosonic and fermionic systems across the transition from integrabil- 
ity to quantum chaos. These systems are integrable in the presence of only nearest-neighbor terms, 
whereas the addition of next-nearest neighbor hopping and interaction may lead to the onset of 
chaos. We show that the strength of the next-nearest neighbor terms required to observe clear sig- 
natures of nonintegrability is inversely proportional to the system size. Interestingly, the transition 
to chaos is also seen to depend on particle statistics, with bosons responding first to the integrability 
breaking terms. In addition, we discuss the use of delocalization measures as main indicators for 
the crossover from integrability to chaos and the consequent viability of quantum thermalization in 
isolated systems. 
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I. INTRODUCTION 



Random matrix theory (RMT) deals with the statisti- 
cal properties of ensembles of matrices composed of ran- 
dom elements. It was originally designed by Wigner in 
his efforts to understand the statistics of energy levels 
of nuclei [1] and was further elaborated by several au- 
thors, notably Mehta p^. RMT received a significant 
boost with the discovery of its connection with classical 
chaos In particular, it was observed that quantum 

systems whose classical analog are chaotic show the same 
fluctuation properties predicted by RMT. 

The application of RMT was soon extended to the de- 
scription of other quantum many-body systems, such as 
atoms, molecules, and quantum dots |7l-[lQ|, and it was 
not restricted to statistics of eigenvalues but accommo- 
dated also the analysis of eigenstates [Tll-[T3| . Important 
developments that led to the broadening of the theory in- 
clude the introduction of ensembles of random matrices 
that take into account the predominance of short range 
interactions in real many-body systems p!3-[l6|, the inti- 
mate connection between quantum trans^rt and spec- 
tral properties of mesoscopic systems [8, \^ [l8| , and the 
relationship between chaos and quantum thermalization 



It has been conjectured that the thermalization of fi- 
nite isolated quantum systems is closely related to the 
onset of chaos and occurs at the level of individual states 
0, [lO, mi 5 which has become known as the eigenstate 
thermalization hypothesis (ETH). Related work was done 
with nuclear shell model calculations and delocalization 
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measures [12|, |26|. More recently, this subject has re- 
ceived renewed attention due to its relevance to ultra- 
cold gas experiments. For example, in a remarkable ex- 
periment by a group at Penn State [l^, it was shown 
that, after being subject to a strong perturbation, a gas 
of bosons trapped in a (quasi-)one-dimensional geometry 
(created by means of a deep two-dimensional optical lat- 
tice) did not relax to the standard prediction of statistical 
mechanics. In contrast to those results, in another exper- 
iment in which a bosonic gas was trapped in a different 
(quasi-) one-dimensional geometry (generated by an atom 
chip), relaxation to a thermal state was inferred to occur 
in a very short time scale [28| . 

Following those experiments, several theoretical works 
have explored the question of thermalization in noninte- 
grable isolated quantum systems after a quantum quench 
in one dimension poMs?} . After numerically exploring the 
nonequilibrium dynamics in finite one dimensional (ID) 
systems, thermalization was observed in some regimes 
[S, [U, 132^ but not in others even though in all 

cases integrability was broken. Several factors may play a 
role in the absence of thermalization in finite ID systems 
after a quench: (i) the proximity to integrable points 
[31I . [32} . (ii) the proximity of the energy of the initial 
nonequilibrium state after the quench to the energy of the 
ground state [3l|, [sH, [s^ , (iii) particle statistics and the 
observable considered (in fermionic systems, the momen- 
tum distribution function may take much longer to relax 
to equilibrium than other observables [sH); and finally 
(iv) quenching the system across a superfluid/metal to 
insulator transition j29|, [30|, [13 • Recent numerical stud- 
ies for bosons and fermions in one dimension have shown 
that there is a direct link between the presence (absence) 
of thermalization and the validity (failure) of the ETH 

In the present work, we provide a detailed description 
of the integrable-chaos transition in the one-dimensional 
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bosonic and fermionic systems studied in Refs. [3l| and 
|32| . These systems are clean and have only two-body in- 
teractions; the transition to chaos is achieved by increas- 
ing the strength of next-nearest-neighbor (NNN) terms 
rather than by adding random parameters to the Hamil- 
tonian. Under certain conditions these systems may also 
be mapped onto Heisenberg spin-1/2 chains. Several pa- 
pers have analyzed spectral statistics of disordered |38l - 
l42j and clean |43l-l46j ID Heisenberg spin-1/2 systems. 
Mostly, they were limited to sizes smaller than consid- 
ered here and, in the case of clean systems, focused on 
properties associated with the energy levels, while here 
eigenvectors are also analyzed. Our goal is to establish 
a direct comparison between indicators of chaoticity and 
the results obtained in Refs. and [32^ for thermaliza- 
tion and the validity of ETH. Our analysis also provides 
a way to quantify points (i) and (ii) in the previous para- 
graph, which can result in the absence of thermalization 
in finite systems. 

Overall, the crossover from integrability to chaos, 
quantified with spectral observables and delocalization 
measures, mirrors various features of the onset of ther- 
malization investigated in Refs. [3l| and [32|, in partic- 
ular, the distinct behavior of observables between sys- 
tems that are close and far from integrability, and be- 
tween eigenstates whose energies are close and far from 
the energy of the ground state. We also find that the 
contrast between bosons and fermions pointed out in 
Ref. [32I is translated here into the requirement of larger 
integrability-breaking terms for the onset of chaos in 
fermionic systems. Larger system sizes also facilitate the 
induction of chaos. In addition, we observe that mea- 
sures of the degree of delocalization of eigenstates become 
smooth functions of energy only in the chaotic regime, a 
behavior that may be used as a signature of chaos. 

The paper is organized as follows. Section HIl describes 
the model Hamiltonians studied and their symmetries. 
Section Uni analyzes the integrable-chaos transition based 
on various quantities. After a brief review of the unfold- 
ing procedure. Sec. III. A focuses on quantities associated 
with the energy levels, such as level spacing distribution 
and level number variance. Section III.B introduces mea- 
sures of state delocalization, namely information entropy 
and inverse participation ratio (IPR), showing results for 
the former in the mean field basis. Results for the inverse 
participation ratio and discussions about representations 
are left to the Appendix. Concluding remarks are pre- 
sented in Sec. |TVl 



II. SYSTEM MODEL 

We consider both scenarios: hardcore bosons and spin- 
less fermions on a periodic one-dimensional lattice in 
the presence of nearest-neighbor (NN) and next-nearest- 
neighbor (NNN) hopping and interaction. The Hamilto- 
nian for bosons Hb and for fermions Hp are respectively 



given by 
Hb = 

L 



[-t + h.c.) +v(^nt- (^nti - \ 

- t' (blbi+2 + h.c)j + V 



,(1) 



and 
Hf 

L 
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Above, L is the size of the chain, hi and h\ {fi and //) 
are bosonic (fermionic) annihilation and creation oper- 
ators on site and = blbi {n{ = fjfi) is the boson 
(fermion) local density operator. Hardcore bosons do not 
occupy the same site, i.e., b]'^ = 6f , so the operators com- 
mute on different sites but can be taken to anti-commute 
on the same site. The NN (NNN) hopping and interac- 
tion strengths are respectively t {t') and V iV'). Here, 
we only study repulsive interactions (V, ^ > 0). We take 
h = 1 and t = y = 1 set the energy scale in the remaining 
of the paper. 

The bosonic (fermionic) Hamiltonian conserves the to- 
tal number of particles (A^/) and is translational in- 
variant, being therefore composed of independent blocks 
each associated with a total momentum k. In the par- 
ticular case of /c = 0, parity is also conserved, and at 
half-filling, particle-hole symmetry is present, that is, the 
bosonic [fermionic] model becomes invariant under the 

transformation flf (^1 + ^i) [Ilf (// + /i)]^ which annihi- 
lates particles from filled sites and creates them in empty 
ones. The latter two symmetries will be avoided here by 
selecting /c 7^ and A'^j = L/3. For even L, we consider 
k = l,2,...,L/2-landfor odd L, k = 1, 2, . . . (L-l)/2. 
The dimension Dk of each symmetry sector studied is 
given in Table HI 

Exact diagonalization is performed to obtain all eigen- 
values and eigenvectors of the systems under investiga- 
tion. When t' = V = models ([1]) and (|2]) are integrable 
and may be mapped onto one another via the Jordan- 
Wigner transformation [[47]. A correspondence with the 
Heisenberg spin-1/2 chain also holds, in which case the 
system may be solved with the Bethe ansatz |49|] . 



III. SIGNATURES OF QUANTUM CHAOS 

The concept of exponential divergence, which is at the 
heart of classical chaos, has no meaning in the quan- 
tum domain. Nevertheless, the correspondence princi- 
ple requires that signatures of classical chaos remain in 
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TABLE I: Dimensions of subspaces 



Bosons 










L = 18 


/c = 1,5,7 


A: = 2,4,8 


/c - 3 


k = 6 




1026 


1035 


1028 


1038 


L = 21 


k = 7 


other /c's 








5538 


5537 






L = 24 


odd /c's 


A: = 2,6, 10 


/c = 4 


/c = 8 




30624 


30664 


30666 


30667 


Fermions 










L = 18 


/c = 1,5,7 


/c = 2,4,8 


/c = 3 


/c = 6 




1035 


1026 


1038 


1028 


L = 21 


k = 7 


other /c's 








5538 


5537 






L = 24 


odd /c's 


/c = 2,6, 10 


A: = 4 


k^S 




30624 


30664 


30667 


30666 



the quantum leveL Different quantities exist to identify 
the crossover from the integrable to the non-integrable 
regime in quantum systems. We consider both spectral 
observables associated with the eigenvalues and quanti- 
ties used to measure the complexity of the eigenvectors. 



A. Spectral observables 

Spectral observables, such as level spacing distribution 
and level number variance are investigated below. They 
are intrinsic indicators of the integrable-chaos transition. 
Their analysis are based on the unfolded spectrum of each 
symmetry sector separately. 



1. Unfolding procedure 

The procedure of unfolding consists of locally rescal- 
ing the energies as follows. The number of levels with 
energy less than or equal to a certain value E is given 
by the cumulative spectral function, also known as the 
staircase function, N{E) = 0(£^ — En), where 6 is 
the unit step function. N{E) has a smooth part Nsm{E), 
which is the cumulative mean level density, and a fluctu- 
ating part Nfi{E), that is, N{E) = Nsm{E) + Nfi{E). 
Unfolding the spectrum corresponds to mapping the en- 
ergies {£^1, £^2, . . . , Ed} onto {ei, 62, . . . e^}, where Cn = 
Nsm{En), SO that the mean level density of the new 
sequence of energies is 1. Different methods are used 
to separate the smooth part from the fluctuating one. 
Statistics that measure long-range correlations are usu- 
ally very sensitive to the adopted unfolding procedure, 
while short-range correlations are less vulnerable (Hof . 
Here, we discard 20% of the energies located at the edges 
of the spectrum, where the fluctuations are large, and 
obtain Nsm{E) by fitting the staircase function with a 
polynomial of degree 15. 



2. Level spacing distribution 

The distribution of spacings s of neighboring energy 
levels [2I, 0, m, [lo| is the most frequently used observ- 
able to study short-range fluctuations in the spectrum. 
Quantum levels of integrable systems are not prohib- 
ited from crossing and the distribution is Poissonian, 
Pp{s) = exp(— s). In non-integrable systems, crossings 
are avoided and the level spacing distribution is given by 
the Wigner-Dyson distribution, as predicted by random 
matrix theory. The form of the Wigner-Dyson distribu- 
tion depends on the symmetry properties of the Hamilto- 
nian. Ensembles of random matrices with time reversal 
invariance, the so-called Gaussian orthogonal ensembles 
(GOEs), lead to Pwd{s) = (7r5/2) exp(-7rsV4). The 
same distribution form is achieved for models ^ and (|2]) 
in the chaotic limit, since they are also time reversal in- 
variant. However, these systems differ from GOEs in the 
sense that they only have two-body interactions and do 
not contain random elements. Contrary to GOEs and to 
two-body random ensembles [l5|, the breaking of sym- 
metries here is not caused by randomness, but instead 
by the addition of frustrating next-nearest-neighbor cou- 
plings. Notice also that the analysis of level statistics in 
these systems is meaningful only in a particular symme- 
try sector; if different subspaces are mixed, level repulsion 
may be missed even if the system is chaotic [46|, [51| . 




FIG. 1: (Color online.) Level spacing distribution for hard- 
core bosons averaged over all /c's in Table [H for L = 24, and 
t^ = V\ For comparison purposes, we also present the Poisson 
and Wigner-Dyson distributions. Bottom right panel: energy 
difference between first excited state Ei and ground state Eo 
in the full spectrum times L, for L = 18 (circles), L = 21 
(squares), and L = 24 (triangles). 

In Figs. (dJ) and (j2j), we show P{s) across the transi- 
tion from integr ability to chaos for bosons and fermions, 
respectively, in the case of L = 24. An average over all 
/c's is performed, but we emphasize that the same be- 
havior is verified also for each /c-sector separately. As 
t\V^ increases and symmetries are broken, level repul- 
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sion becomes evident, the peak position of the distribu- 
tion shifts to the right, and the tail of the distribution 
changes from exponential to Gaussian. Excellent agree- 
ment with the Wigner-Dyson distribution is seen already 
for t' = V' > 0.12. The bottom right panels in Figs. ([I]) 
and ([2]) give the energy difference between first excited 
state and ground state times L as a function of t\V\ 
One can see there that the product is size independent 
emphasizing that the ground state of the systems con- 
sidered here is gapless in the thermodynamic limit, as 
expected from the phase diagrams presented in Ref. [52| . 
Notice that the particular case of the fermions exhibits 
an even-odd finite-size effect that becomes irrelevant in 
the thermodynamic limit. 



0.15 




FIG. 2: (Color online.) As in Fig.[T]but for spinless fermions. 

To better quantify the integrable-chaos transition, we 
show in Fig. [3] the level spacing indicator a, defined as 
follows: 



(3) 



where the sums runs over the whole spectrum. We should 
stress that this is a discrete rather than integral sum, be- 
cause P{s) as computed by us is a discrete quantity. For 
a chaotic system a ^ 0. The indicator a is comparable 
to the quantity r] introduced in Ref. 



22|. 



As seen from the bottom panels and insets in Fig.[3l the 
values of t', leading to the transition to chaos decrease 
with the size of the system, suggesting that the onset 
of chaos in the thermodynamic limit might be achieved 
with an infinitesimally small integrability breaking term, 
although the existence of a saturation value cannot be 
discarded [45|. A conclusive statement would require 
even larger systems or a theory for the behavior of sys- 
tems approaching infinite sizes. Interestingly, a decays 
faster for bosons, which indicates that the crossover to 
the chaotic behavior may depend on particle statistics. 
This contrasts studies of the ground state properties of 
many-body systems with two-body interactions, where 
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FIG. 3: (Color online.) Average a over all /c's; t' = V' . Top 
panel: difference a (fermions) - a (bosons). Left bottom panel: 
bosons; right bottom panel: fermions. Semi- logarithmic plot 
in main panels and linear plot in insets. Circles: L — 18, 
squares: L = 21, triangles: L = 24. 



such differences were not found [53|. The top panel in 
Fig. [3] shows the difference between the value of a for 
fermions and bosons. In general, it diminishes with in- 
creasing the size of the chain, and the point at which the 
difference attains its maximum value moves toward lower 
values of t' . 




FIG. 4: (Color online) Peak position of the level spacing 
distribution averaged over all /c's; t' ^V'. Top panel: dif- 
ference between the peak position of bosons and fermions, 
peak(bosons) - peak (fermions). Left bottom panel: bosons; 
right bottom panel: fermions. Semi-logarithmic plot in main 
panels and linear plot in insets. Circles: L — 18, squares: 
L = 21, triangles: L — 24. Dashed line indicates the peak 
position of Pwd{s). 

The findings in Fig. [3] are reinforced in Fig.lH where we 
show the approach of the peak position of the level spac- 
ing distribution to the peak position of Pwd{s) as t' ^ V 
increases. The transition is faster for larger chains and 
once again depends on particle statistics, with bosons re- 
sponding first to the breaking of symmetries. We should 
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add that equivalent results are obtained by fitting P{s) 
with the Brody distribution [l5|, 



Pb{s) = exp , b ■ 



^+1 



and analyzing how the increase in t\ V changes /3 from 
(in the integrable region) to 1 (in the chaotic limit). 

The spectral properties of models ([1]) and (|2]) discussed 
here are somehow mirrored by their dynamical behavior, 
which were studied respectively in Refs. [31] and (32j . 
(Notice that the analysis of the integrable-chaos transi- 
tion performed here uses the same values of t' con- 
sidered in those works.) The smooth approach to inte- 
grability, as shown in the Figs. 1-4 above, is followed by 
the breakdown of thermalization observed in Refs. |31j 
and [32|. However, an important difference between our 
results here and the results in Refs. [3l| and [32| is that 
in the latter works it was not clear that the values of 
t', V required for the system to thermalize would reduce 
with increasing system size, whereas this is the case here 
for obtaining a Wigner-Dyson distribution of the level 
spacings. 

Another interesting feature found in Ref. [32[ is that 
in the context of quenched dynamics there are differ- 
ences associated with the particle statistics. In partic- 
ular, it was seen that some observables such as the mo- 
mentum distribution function [n(/c)] in fermionic systems 
may take longer time to relax to equilibrium than their 
bosonic counterparts. (A related effect in which a quasi- 
steady regime occurs for n(/c) before full relaxation has 
been suggested for higher dimensional fermionic systems 
[53-[56|.) It was also shown in Ref. [32| that the differ- 
ence between the eigenstate expectation values of n(/c) 
for eigenstates of the fermionic Hamiltonian with close 
energies suffer from particularly large finite size effects 
when compared to other observables such as the density- 
density structure factor [A^(/c)] and when compared to 
n(/c) and A/'(/c) for hardcore bosons. Interestingly, here 
we find that for these finite size systems the measures 
of chaoticity also exhibit differences between hardcore 
bosons and fermions, where the former ones respond first 
to the integr ability breaking terms. Further studies will 
be required to explore the relation between the latter 
finding and the onset of ETH for different observables 
and different particle statistics in ID systems. 



S. Level number variance 

Other quantities sensitive to spectral fluctuations in- 
clude measures of long-range correlations, such as spec- 
tral rigidity and level number variance [8]. Both are 
closely related and measure the deviation of the stair- 
case function from the best fit straight line. Here, we 
show results for the level number variance, 5]^(/), defined 
as 



where N{l,e) gives the number of states in the interval 
[e, e+/] and (.) represents the average over different initial 
values of e. For a Poisson distribution, S^(/) = /, while 
for GOEs in the limit of large /, E2(/) = 2[ln(27r/) + 
7 + 1 — 7r^/8]/7r^, where 7 is the Euler constant. Level 
repulsion leads to rather rigid spectra and fluctuations 
become much less significant than in the random energy 
sequences of regular systems. 




E'il)^{iNil,er)-{Nil,e)y 



(4) 



FIG. 5: (Golor online.) Level number variance averaged over 
all /c's. In each panel, solid lines from top to bottom: t' — 
V = 0,0.02,0.04,0.08,0.16,0.32,0.64. Dashed line: GOE, 
dotted-dashed line: Poisson. 

As expected, the level number variance shown in Fig. [5] 
approaches the GOE curve as the strength of NNN in- 
teractions increases. The proximity to the GOE result 
also improves as the system size increases. However, 
deviations from the GOE curve are verified for values 
of t' where the level spacing distribution is already 
very close to a Wigner-Dyson, especially in the case of 
fermions (right panels). In fact, the distinct behavior 
associated with particle statistics becomes yet more ev- 
ident when studying The level number variance 
for the bosonic system with L = 24, for example, coin- 
cides with the GOE result for a large range of values of 
/ already when t' = 0.32 and 0.64, whereas the same 
is not verified for fermions. This further supports the 
view that particle statistics may play an important role 
in the relaxation dynamics and thermalization of finite 
isolated quantum systems. On the other hand, the size 
dependence of the results is an indicator that in the ther- 
modynamic limit the differences between the quantities 
discussed in this paper may become negligible when com- 
paring bosons and fermions, a conjecture that deserves 
further investigation together with its implications for 
the dynamics and thermalization of those systems. 



B. Delocalization measures 

Contrary to spectral observables, quantities used to 
measure the complexity of eigenvectors, as delocaliza- 
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tion measures [ll|, [l2| , are not intrinsic indicators of the 
integrable-chaos transition since they depend on the ba- 
sis in which the computations are performed. The choice 
of basis is usuahy physicahy motivated. The mean-field 
basis is the most appropriate representation to separate 
global from local properties, and therefore capture the 
transition from regular to chaotic behavior [12|. Here, 
this basis corresponds to the eigenstates of the integrable 
Hamiltonian (t^, = 0). Other representations may also 
provide relevant information, such as the site basis, which 
is meaningful in studies of spatial localization, and the 
momentum basis, which can be used to study /c-space 
localization (see the Appendix for further discussions). 

The degree of complexity of individual eigenvectors 
may be measured, for example, with the information 
(Shannon) entropy S or the inverse participation ratio 
(IPR). The latter is also sometimes referred to as number 
of principal components. For an eigenstate ipj written in 



the basis vectors (j)k as ipj 
respectively given by 



D 

■El 

k=l 



In 



S and IPR are 



(5) 



and 



IPR, 



E 



(6) 



The above quantities measure the number of basis vectors 
that contribute to each eigenstate, that is, how much 
delocalized each state is in the chosen basis. 

For the GOE, the amplitudes Cj are independent ran- 
dom variables and all eigenstates are completely delocal- 
ized. Complete delocalization does not imply, however, 
that S = lni:>. For a GOE, the weights |c^p fluctu- 
ate around 1/D and the average over the ensemble is 
Sgoe = ln(0.48I)) + (9(1/1^) MM- 

Figures [6] and [7] show the Shannon entropy in the mean 
field basis Sm/ vs the effective temperature for bosons 
and fermions, respectively. The effective temperature, 
Tj of an eigenstate tjjj with energy Ej is defined as 



where 



Z = Tr {e-^/^^ } . 



(7) 



(8) 



Above, H is Hamiltonian ([T]) or (j2j), Z is the partition 
function with the Boltzmann constant /c^ = 1, and the 
trace is performed over the full spectrum as in Refs. (sij 
and [32I (see the Appendix for a comparison with effec- 
tive temperatures obtained by tracing over exclusively 
the sector k = 2). The figures include results only for 
Tj < 10; for high Ej^ the temperatures eventually become 
negative. By plotting the Shannon entropy as a function 
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FIG. 6: (Color online.) Shannon entropy in the mean field 
basis vs effective temperature for bosons, L = 24, k = 2, 
and = V\ The dashed line gives the GOE averaged value 
Sgoe - ln(0.48L>). 
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FIG. 7: (Color online.) As in Fig. [6] for fermions. 



of the effective temperature, we allow for a direct com- 
parison of our results here and the results presented in 
Refs. [lij and m. 

As seen in Figs. [6] and [71 the mixing of basis vectors, 
and therefore the complexity of the states, increases with 
t', y , but it is only for Tj > 2 that the eigenstates of our 
systems approach the GOE result. Similarly, in plots of 
Sm/ vs energy (see Figs. [TT] and [12] in the Appendix), it is 
only away from the borders of the spectrum that Smf 
Sgoe; in the borders, the states are more localized and 
therefore less ergodic. This feature is typical of systems 
with a finite range of interactions, such as models ([1]) 
and ([2]) and also banded, embedded random matrices, 
and two-body random ensembles [l3|, HH, 113 • 

The analysis of the structure of the eigenstates hints 
on what to expect for the dynamics of the system. In 
the context of relaxation dynamics, not only the den- 
sity of complex states participating in the dynamics is 
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relevant, but also how similar these states are. In close 
connection, but from a static perspective, it is the on- 
set of chaos that guarantees the uniformization of the 
eigenstates. According to Percival's conjecture [58|, the 
complexity of chaotic wave functions adjacent in energy is 
very similar, they essentially show the same information 
entropy. A further extension of this idea is Berry's con- 
jecture (Hil, which assumes that energy eigenfunctions in 
a time-reversal invariant and ergodic system is a super- 
position of random plane waves. The Eigenstate Ther- 
malization Hypothesis (ETH) [20| can be related to the 
vahdity of Berry's conjecture. ETH states that thermal- 
ization of an isolated quantum system occurs when each 
eigenstate already exhibits a thermal value for the ob- 
servables, that is, the eigenstate expectation values do 
not fluctuate between eigenstates close in energy. 

In Figs. [6] and [7] (see also Figs. [TT]fT6l in the Appendix), 
the structure of the eigenstates close in energy reveals 
fluctuations throughout the spectrum as we approach in- 
tegr ability, whereas in the chaotic regime, fluctuations 
are mostly restricted to the edges of the spectrum, with 
Smf being a smooth function of energy (or temperature) 
away from the borders. A related result was seen in 
Ref. [4l| , where a clear relationship between an entangle- 
ment measure and a delocalization measure for a clean 
Heisenberg model appeared only in the chaotic limit (two 
dimensions), being absent in its integrable counterpart 
(one dimension). Fluctuations imply that eigenstates 
very close in energy have different degrees of complexity 
and localization properties; these states may therefore 
not entirely comply with the ETH. This explains the ab- 
sence of thermalization in the integrable limit. In fact, 
a similar conclusion was achieved in Refs. [3l| and [32| 
based on plots for the momentum distribution function 
vs energy (Figs. 4 and 7, respectively). There, it was 
observed that fluctuations between expectation values of 
states close in energy increase toward integrability. 



T ' 1 ' 1 1 I 1 ' 1 ' r 




E E 

FIG. 8: (Color online.) Shannon entropy in the mean field 
basis vs energy when the system is close to integrability (top 
panels) and in the chaotic limit (bottom panels), k — 2, and 
= V' . Panels on the left: bosons; panels on the right: 
fermions. Curves from bottom to top: L — 18, 21, 24. 

Two additional points need to be made here. First, 



even in the chaotic regime one sees that states in the 
edge of the spectrum remain "localized" in the mean-field 
basis. This is accompanied by a failure of ETH in that 
regime [111, [32| , and hence thermalization is not expected 
to occur when the energy of the time-evolving state af- 
ter the quench is close to the ground-state energy (low 
effective temperatures). This is an important feature of 
isolated quantum systems that will need to be considered 
with more care when dealing with ultracold gases ex- 
periments. Second, in comparing bosons and fermions, 
larger fluctuations are verified for the latter, offering a 
further justification for the deviations between statisti- 
cal mechanics predictions for observables after relaxation 
and the exact time-averaged result of the quantum evo- 
lution observed in finite fermionic systems Note, 
however, that fluctuations appear to decrease with sys- 
tem size, as shown in Fig. [8] (specially noticeable in the 
bottom right panel). This may be a simple reflection of 
better statistics, but may suggest also that in the thermo- 
dynamic limit some of the differences between fermions 
and bosons may eventually disappear. 



IV. CONCLUSIONS 

We have presented a detailed analysis of the transi- 
tion from integrability to quantum chaos for gapless one- 
dimensional systems of interacting spinless fermions and 
hardcore bosons. Here, the onset of chaos was dictated by 
the enhancement of next-nearest-neighbor hopping and 
interactions. 

Our comparisons for different system sizes suggested 
that in the thermodynamic limit an infinitesimal inte- 
grability breaking term suffices for the onset of chaos, 
although further studies are necessary for settling this is- 
sue. Also, this may not warrant that thermalization will 
occur for infinitesimal integrability breaking terms since, 
at least for our finite systems, we could not establish a 
one to one correspondence between the two effects. 

We have found differences in behaviors associated with 
particle statistics. The transition to chaos in fermionic 
systems, as measured by level spacing indicator, peak po- 
sition of level spacing distribution, and level number vari- 
ance, required integr ability-breaking terms larger than in 
the bosonic case. With respect to delocalization mea- 
sures, larger fluctuations were also verified for fermions. 

We studied wave function complexity using different 
delocalization measures and choices of underlying ba- 
sis. Our results have shown that the similar structure of 
eigenstates close in energy is a primary feature of chaotic 
systems. This finding reinforces the proposal to elevate 
Berry's conjecture to the status of the best definition 
of quantum chaos [13] and suggests that the onset of a 
smooth dependence of delocalization measures with en- 
ergy be used as an indicator of quantum chaos and a 
condition for quantum thermalization. 

Finally, we have shown that even when the systems 
are chaotic in terms of the level spacing distribution and 
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the level number variance, there are still regions in the 
edges of the spectrum in which the states are less de- 
localized and their structures less similar. As shown in 
Refs. [sil and [32[, those states do not satisfy ETH and 
hence, whenever one performs a quench in a system so 
that the energy of the time evolving state is close to the 
ground state energy [or in other words, when the effective 
temperature of the system as defined by Eq. ([7]) is very 
low] , relaxation of observables to the thermal distribution 
prediction is not expected. 

The analysis and findings described here are intimately 
refiected b y t he studies of thermalization pursued in 
Refs. [ill, jsij . and [32| and provide strong support to 
those works. Moreover, the lattices we considered may 
also be mapped onto other one-dimensional systems, such 
as spin- 1/2 chains, which indicates the broad range of ap- 
phcability of our results for gapless systems. 



Acknowledgments 

L.F.S. thanks support from the Research Corporation. 
M.R. was supported by the US Office of Naval Research 
and by Georgetown University. 



10 







- Bosons 


/ / 1 ' 











- Fermions 


' 11' 







-20 



-15 -10 -5 
Energy 



-20 



-15 -10 -5 
Energy 



FIG. 9: (Color online.) Taii vs energy, where Taii stands for 
the effective temperature computed considering the eigenval- 
ues of all symmetry sectors; L = 24. Curves from bottom to 
top: t' = V' = 0.00, 0.02, 0.03, 0.04, 0.06, 0.08, 0.12, 0.16, 
0.24, 0.32, 0.48, 0.64, 0.96, 1.28. 



in the scale of Fig. [10] due to signiffcant discrepancies 
between Taii and Tk=2' 

More generally, for larger system sizes, the tempera- 
ture is expected to be computed by means of quantum 
Monte-Carlo simulations or other better scaling numeri- 
cal approaches. This means that in general all sectors will 



Appendix A: Complexity of the wave functions 

We provide here further illustrations for the complex- 
ity increase of the wave functions with the onset of chaos 
for models ([T]) and (j2j). This is based on the computation 
of Shannon entropy, Eq. ([5]), and the inverse participa- 
tion ratio, Eq. (|6]). We compare results in both repre- 
sentations, mean-ffeld- and /c-space-bases. Overall, the 
approach to chaos is followed by the reduction in ffuctua- 
tions in the results for S and IPR close in energy, with the 
decrease in ffuctuations being slower for fermions than for 
bosons. 



1. Effective temperature 

Figures [6] and [3 gave the entropy in the mean-ffeld basis 
vs the effective temperature. There, each temperature 
Tj , for an eigenstate of energy Ej , was obtained by means 
of Eq. ([7j) and performing the trace over the full spectrum 
(let us call it Taii here). In Fig. [9l we compare Taii with 
the eigenstate energies. 

One may also wonder what would happen if one uses 
only the spectrum of the k = 2 sector to perform the trace 
and hence to compute the temperature T/e=2- Actually, 
for the system sizes employed here, the values obtained 
in this latter way do not differ much from temperatures 
calculated considering the energies of all /c-sectors [g^. 
Figure [10] shows that the largest disagreements between 
Tall and Tk^2 occur at low energies; but even then, they 
are usually not higher than 5%. Exceptions are the ffrst 
two or three lowest temperatures, which do not appear 
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FIG. 10: (Color online.) Temperature difference vs Taii- 
Tdiff = 100\Taii -Tk^2\/Taii, where Tfc=2 is computed consid- 
ering only the eigenvalues from the k = 2 sector. Left panels: 
bosons; right panels: fermions; t' = V' . Negligible differences 
are seen between the curves for L = 24 [light gray (green - 
online)] and L = 21 [dark gray (red - online)] when Taii > 2, 
while the curve for L = 18 (black) saturates at a higher level 
(especially noticeable for fermions). 
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be considered when computing T. Our results here show 
that the differences with considering specific momentum 
sectors are small and decreasing with the system size. 



2. Results in the mean-field basis 

Figures [TT] and [12] show the mean-field Shannon en- 
tropy vs energy for all the eigenstates of the k = 2 sector. 
Notice that the whole spectrum for the k = 2 sector is 
presented and not only the energies leading to T < 10 
as in Figs. [6] and [71 The typical behavior of banded ma- 
trices is observed: larger delocalization appearing away 
from the edges of the spectrum, although not as large as 
the GOE result Sqoe = MOASD)-h 0{1/D), and lower 
complexity at the edges [isl, UH, |53 • 
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FIG. 11: (Color online.) Shannon entropy in the mean field 
basis vs energy for bosons, L = 24, /c = 2, and = V' . The 
dashed line gives the GOE averaged value Sgoe ^ ln(0.48D). 
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FIG. 12: (Color online.) Same as in Fig. [TT] for fermions. 

A similar behavior is seen in the plots of the in- 
verse participation ratio in the mean-field basis vs en- 



ergy (Figs. [13] and [T4|). The IPR values increase sig- 
nificantly with y^ but do not reach the GOE result 
IPR ={D^ 2)/3 [m, [H. IPR gives essentially the same 
information as S, although the first shows larger fiuctua- 
tions. 
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FIG. 13: (Color online.) Inverse participation ratio in the 
mean field basis vs energy for bosons, L = 24, /c = 2, and 
t' = V'. The GOE result IPRgoe - i:>/3 is beyond the 
chosen scale. 
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FIG. 14: (Color online.) Same as in Fig. [13] for fermions. 



3. Results in the /c— basis 

Identifying the mean-field basis may not always be a 
simple task. For example, some ID models may have 
more than one integrable point. It may also happen that 
one is so far from any integrable point that there is no rea- 
son to believe that such a point has any relevance for the 
chosen system. The latter case may be particularly ap- 
plicable to higher-dimensional systems where integrable 
points are, in general, the noninteracting limit or other 
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trivial limit. Working on the mean- field basis also adds 
an extra step in the computations since the diagonal- 
ization of the system is usually not performed in that 
basis, i.e., one needs to perform a change in basis when 
computing S and IPR in the mean-field basis. This ex- 
tra computation step may become very demanding when 
dealing with large systems. In addition, depending on 
the studies being performed, it may be of interest to an- 
alyze the structure of the eigenvectors in another basis. 
The problem of spatial localization, for example, calls for 
the use of the site-basis. 
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FIG. 15: (Color online.) Shannon entropy in the /c-basis vs 
energy for bosons, L = 24, k = 2, and = V' . The dashed 
line gives the GOE averaged value Sgoe ^ ln(0.48D). 

Motivated by the discussion above, we include here 
the results for S in the /c-basis in which we diagonalize 
our Hamiltonians. Those are shown in Figs. [15] and 
[T6l Large entropy values are now simply related to 
high delocalization with respect to the /c-basis and 
have nothing to do with the onset of chaos. They are 
found in both integrable and non-integrable regimes 



and may even surpass Sgoe- Other differences between 
the mean-field and k bases include: (i) the localization 
increase expected for both edges of the spectrum in 
banded matrices is not so evident in the /c-basis, some 
high energy states remaining as delocalized as the 
central states; and (ii) the distinct degree of fiuctuations 
between bosons and fermions, even though still higher 
for fermions, is not so visible anymore. In spite of these 
deviations, the /c-basis may still be used as a signature of 
the integr able-chaos transition. The reason being that, 
just as in the mean field basis, the dependence of S^ 
with energy becomes smoother only in the chaotic limit. 
Therefore, since reduction in fiuctuations in S and IPR 
for states close in energy has been pointed as a main 
cause for the validity of the ETH, the /c-basis may still 
be used to determine where the onset of thermalization 
is expected. 
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FIG. 16: (Color online.) Same as in Fig. [15] for fermions. 
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